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Abstract 

Measurements are presented of meson and A baryon production in deep-inelastic 
positron-proton scattering (DIS) in the kinematic range 10 < < 70 GeV^ and 
10~^ < a; < 10~^. The measurements, obtained using the HI detector at the HERA 
colhder, are discussed in the hght of possible mechanisms for increased strangeness 
production at low Bjorken-x. Comparisons of the xp spectra, where Xi;' is the frac- 
tional longitudinal momentum in the hadronic centre-of-mass frame, with results 
from electron-positron annihilation are made. The xp spectra and the "seagull" 
plot are compared with previous DIS results. The mean and A multiplicities 
are studied as a function of the centre-of-mass energy W and are observed to be 
consistent with a logarithmic increase with W when compared with previous mea- 
surements. A comparison of the levels of strangeness production in diffr active and 
non-diffractive DIS is made. An upper limit of 0.9 nb, at the 95% confidence level, 
is placed on the cross-section for QCD instanton induced events. 
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1 Introduction 



This paper describes studies of the flavour structure of the hadronic final state in deep-inelastic, 
neutral current, positron-proton scattering (DIS) performed using the HI detector at the HERA 
collider at DESY. The scattering process is described in terms of the variables = —q'^, where 
q is the four-momentum of the exchanged boson, x = Q'^/2p ■ q and y = p • q/p • k where p 
and k are, respectively, the four-momenta of the initial state proton and electron. Results from 
HERA have shown that the proton structure function F2{x,Q'^) measured in DIS rises rapidly 
as the Bjorken variable x decreases below about x = 10~^ [||, ^. Here it is examined if this rise 
is accompanied by a change in the flavour composition of the hadronic final state, in particular 
the strangeness production rate, associated with either the hard interaction or the subsequent 
hadronisation process. 

Typically in electron-positron annihilation and DIS experiments at moderate or large x the 
production ratios of u, d, and s quarks in the soft hadronisation process have been found to 
be 1 : 1 : 0.3 |^, ^, that is = 0.3, where is the strangeness suppression factor, the ratio 
of s quark to n or to c? quark production necessary to describe the data using the Lund string 
hadronisation model Q. The partial suppression of s quark production in the hadronisation is 
a consequence of the relative masses of the u, d and s quarks; the mass of the c quark is such 
that c quarks are not produced in the hadronisation process. Most recent DIS results [^, |^, ^ 
favour a suppression of strangeness production stronger than that above. As « 0.2, as does one 



older measurement , although there is a recent result which prefers the higher value of A^ [ IC ] . 



Recent electron-positron annihilation measurements |Tl[ also tend to favour A^ ~ 0.2, but again 



the experimental situation is not clear [12| 



In hadron-hadron collisions the mean value of the strangeness suppression factor is about 
0.3 but there is evidence for a dependence on the centre-of-mass energy, the proportion of 
strangeness produced rising with energy. There is also some evidence for a dependence on the 



region of phase space investigated [|13|. 



HERA offers the possibility of studying strangeness production in situations similar to those 
prevailing in all the above processes. Hadron production in DIS in the phase space region 
associated with the struck quark is expected to resemble that in electron-positron annihilation. 



This is confirmed by recent measurements |TJ, 15]. The rapidity region between the struck 
quark and the proton remnant in DIS is likely to be similar to the central plateau in hadron- 
hadron collisions. A recent comparison of hadron production in photoproduction and deep- 
inelastic electron-proton scattering supports this possibility The proton remnant is likely 
to be similar to the hadron remnants in hadron-hadron scattering. The phase space region in 
which strange particles may be identified using the current HI apparatus corresponds to the 
struck quark and central plateau regions. The broad agreement of the A^ results obtained in 
electron-positron annihilation with the measurements in the struck quark and central plateau 
regions in DIS and the central plateau in hadron-hadron collisions suggests that the colour 
forces responsible for hadron production in these regions are similar. It is of interest to examine 
whether this similarity persists at the lowest values of x accessible at HERA. 

A possible source of strangeness production levels in DIS at HERA above those expected 
from previous measurements is QCD instantons [l^]. Instantons are expected features of both 
QCD and the Glashow-Weinberg-Salam (GWS) model. While the cross-section for instanton 
induced processes is subject to large uncertainties, the properties of the final states of instanton 
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induced events may be calculated. These include exciting signatures such as baryon plus lepton 
number non-conservation in the case of the GWS model and, analogously, non-conservation of 
chirality in the case of QCD. The latter signature is difficult to observe experimentally, but a 
further property of QCD instanton induced events, the democratic production of all kinemat- 
ically allowed flavours, would lead to a measureable increase in the quantity of strangeness in 
the final state. The kinematics of instanton induced events is such that this excess strangeness 
would be produced largely in the regions of the HI detector close to the proton remnant. 

A further possible mechanism for the production of strangeness at levels above expectations 
based on previous measurements may come into play at low x. Although the behaviour of 
F2{x, Q"^) in this region can be accurately described within the framework of perturbative QCD 
using the Dokshitzer, Gribov, Lipatov, Altarelli and Parisi (DGLAP) evolution equations 
it has been argued that the sharp rise in F2{x^Q'^) with decreasing x is an indication that 
new QCD dynamics, described by the Balitsky, Fadin, Kuraev and Lipatov (BFKL) evolution 
equation is starting to play a role at the lowest values of x accessible to HERA [20|. If 



the BFKL description is appropriate it is expected that the hadronic final state will also exhibit 
features associated with BFKL evolution, in particular the production of an increased number of 
gluons with large momentum transverse to the proton in the region around the proton remnant 



direction |21|. There are some indications that this is indeed observed |22|. This would then 
imply an increased level of particle, including strange particle, production in this region. Indeed, 
the proportion of strange particles observed might also increase as there is some evidence that 
events containing energetic gluons in electron-positron annihilation contain a higher proportion 
of mesons than is observed in two-jet events [23]. 



Studies at HERA have revealed that about 10% of DIS events are diffractive in nature [24, 25]. 



These events exhibit a large region in rapidity, close to the proton remnant, in which there is 
no significant hadronic energy. It has been shown that a large proportion of them may be 
described in terms of deep-inelastic positron-pomeron scattering. Furthermore, measurements 
are in agreement with Monte Carlo simulations in which the pomeron is assumed to be an 
object with partonic sub-structure and suggest that the pomeron is composed predominantly of 



gluons [26 



Describing the kinematics of diffractive deep-inelastic scattering (DDIS) using the positron- 
pomeron scattering picture requires the introduction of two further variables: t, the square of 
the four-momentum of the pomeron and xjp, the proportion of the proton's four-momentum 
carried by the pomeron. The latter may be determined using the relationship 



where is the squared invariant mass of the photon-pomeron system and W"^ is the squared 
invariant mass of the complete hadronic final state. The approximation is good as long as both 
the squared proton mass and \ t\ are much smaller than either or W"^^ which is the case in the 
kinematic region considered here. In this picture of DDIS, as the pomeron has vacuum quantum 
numbers, the production of flavour is expected to be democratic up to effects arising from the 
quark masses. 

An alternative model for DDIS is based on the idea that the underlying process is electron- 
gluon scattering via boson-gluon fusion (BGF) followed by the randomisation of the colours 
of the produced quark and anti-quark |2M|. The probability of this producing a colour singlet 
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configuration is g. In this case there is no colour connection between the proton remnant and the 
quark and anti-quark, no hadrons are produced in this region, and the above difFractive signature 
is observed. In this model quark flavour production is again democratic up to kinematic effects. 
It is of interest to examine whether the quantity of strangeness produced in DDIS events is 
compatible with these models of diffraction. 

Also of interest is the possibility that particles containing strange quarks produced in the 
hard interaction may be identified. Possible interaction mechanisms are: BGF resulting in 
the production of an ss pair, cc production via BGF followed by decay to s quarks, and the 
direct interaction of the virtual boson with a strange sea quark. This may ultimately allow 
measurement of the gluon and strange sea contributions to proton structure. 

The above considerations provide the motivation for the work reported here. Strangeness 
production is examined by measuring the numbers and proportions of Kg mesons and A and A 
baryons present in the final state of DIS interactions. In the following, is taken to refer to 
both and and A to both A and A unless explicitly stated to the contrary. 

2 The Detector 

A detailed description of the HI detector can be found in [^]. The following briefly describes 
the components most important to this analysis. The coordinate system used is defined to 
have its origin at the nominal interaction point. The z axis is chosen to be along the proton 
beam direction, the y axis vertical and the x axis completes the right-handed coordinate system. 
Where polar coordinates are used, 9 is the angle with respect to the z axis and (p the azimuthal 
angle with respect to the x axis. Forward and backward refer to the positive and negative z 
directions respectively, hence positive pseudo-rapidities rj = — lntan(0/2) refer to objects in the 
forward hemisphere. 

Surrounding the positron-proton interaction region are detectors designed to measure the 
tracks followed by charged particles. These are the central track detector (CTD) which covers 
the polar angular range 15° < 9 < 165° (2.0 > rj > —2.0) and the forward track detector 
(FTD) covering the range 7° < 9 < 25° (2.8 > rj > 1.5). The CTD consists of concentric 
inner and outer jet drift chambers. Located at the inner radii of these chambers are multi- 
wire proportional chambers (MWPCs), used for triggering purposes, and further drift chambers 
designed to measure accurately the z position of charged particle tracks. A uniform magnetic 
field of 1.15 T parallel to the beam axis is maintained throughout the region containing the 
CTD and FTD by a superconducting solenoid (see below). This system enables the transverse 
momentum of a charged particle to be measured with a resolution of Cpj, /PT ^ O.OOQpT e 0.015 
in the CTD {pr in GeV/c). 

The liquid argon (LAr) calorimeter extends over the polar angular range 4° < 9 < 154° 
(3.4 > rj > —1.5) and covers the entire azimuth. It consists of an electromagnetic section 
with lead absorbers, corresponding to a depth of between 20 and 30 radiation lengths, and a 
hadronic section with steel absorbers. The total depth of the calorimeter varies between 4.5 and 
8 hadronic interaction lengths. 

The backward electromagnetic calorimeter |3^] (BEMC), which is 21.7 radiation lengths 
deep, covers the region 151° < 9 < 176° (—1.4 > r] > —3.4). A resolution of 0.10/^/E with 
E in GeV has been achieved for positrons. The BEMC energy scale is known to an accuracy 
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of 1.7%. Immediately in front of the BEMC is the backward proportional chamber (BPC) 
which determines the point at which the scattered positron enters the BEMC to an accuracy of 
1.5mm over the polar angle range 155° < 6 < 174° (—1.5 > r] > —2.9). In conjunction with a 
measurement of the z coordinate of the primary interaction this leads to an accuracy of 1 mrad 
for the measurement of the polar angle of the scattered positron. 

Behind the BEMC is a time-of-flight system (ToF) with a time resolution of about 1 ns. This 
enables the rejection of background events caused by the interaction of protons in the material 
before the HI detector, as the particles from these events arrive at the ToF system earlier than 
those from interactions at the nominal vertex. 

The LAr calorimeter is surrounded by a superconducting solenoid which provides the mag- 
netic field in the region containing the track detectors. 

Behind the LAr calorimeter and the BEMC is an iron/streamer tube calorimeter which serves 
as the return yoke for the magnetic field and detects hadronic energy leaking out of the inner 
calorimeters. The iron/streamer tube system also serves as a muon detector. 

Surrounding the beam-pipe in the forward direction is a copper/silicon calorimeter, the Plug. 
This device covers the range 0.7° < 6 < 3.5° (5.1 > r] > 3.5). 

Outside the iron in the forward direction is a muon spectrometer which consists of a toroid 
magnet surrounding the beam-pipe sandwiched between two layers of drift chambers. In addition 
to its primary function of detecting forward-going muons this device detects tracks from the 
secondary interactions of particles in the range 0.2° < 6 < 0.8° (6.6 > rj > 5.0). These occur in 
the collimators designed to shield the detector from synchrotron radiation, the forward beam- 
pipe and the adjoining material. 

The luminosity measurement system consists of electromagnetic calorimeters placed at z = 
—33 m and z = —103 m [^]. 



3 Monte Carlo Models 



Calculations of the expected rate of strangeness production in DIS are available in the form of 
Monte Carlo event generators. These are used here, together with a Monte Carlo simulation of 
the HI detector, both to correct the measurements for the finite acceptance and resolution of 
the detector and to make comparisons with theoretical predictions. The standard DIS Monte 
Carlo generators used may be split into two classes. The first class consists of those expected to 
describe the hadronic final state in non-diffractive DIS (NDDIS), but which contain no explicit 
modelling of diffractive DIS. The second class is designed explicitly to model DDIS events. Here 
only one such generator is considered. 

The NDDIS generators are LEPTO 6.1 H, ARIADNE 4.03 |§ and HERWIG 5.8 |3|]. 
These generators simulate positron-proton collisions using the leading order electroweak matrix 
element for the interaction of a virtual boson with a parton, including exact leading order 



QCD corrections. The parton distributions used are the MRS(H) set of the Durham group |36] 
evolved to the appropriate using the DGLAP formalism. These parton distributions were 
determined using most relevant available data, including the early HERA measurements of the 
proton structure function F2{x,Q'^) 38]. The subsequent calculation differs in the three 
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Parameter 


Description 


Default value 


DELPHI value 


E665 value 


Parj(2) 


A. 


0.3 


0.23 


0.2 


Parj(ll) 


Prob. meson has spin 1 (u, d) 


0.5 


0.365 


0.5 


Parj(12) 


Prob. meson has spin 1 (s) 


0.6 


0.410 


0.6 



Table 1: Values of some JETSET hadronisation parameters associated with the production of 
n, d and s quarks required to give a good description of PETRA and PEP data, DELPHI and 
E665 data. 



generators. LEPTO simulates higher order perturbative QCD effects using independent parton 
showers in the leading logarithm approximation to generate radiation from the partons entering 
and leaving the hard sub-process. The observed hadrons are then produced using the Lund 
string hadronisation model |^], as implemented in JETSET 7.4 [^]. Three sets of fragmentation 
parameters are used in JETSET in this study: the default values which give a good description 
of PETRA and PEP electron-positron annihilation data, those providing the best description 



of DELPHI electron-positron annihilation data |11| and those which best describe E665 muon 



nucleon scattering data ||6|. These parameters are listed in table |l|. The term CDM is used to refer 
to a calculation in which higher order perturbative QCD effects are simulated using ARIADNE, 
an implementation of the Colour Dipole Model . In this model gluon radiation takes place 
initially from the dipole formed by the colour charges separated by the hard scattering process 
and subsequently from the dipoles formed by the radiated gluons. Again hadronisation is done 
using JETSET with the above parameter sets. HERWIG simulates parton showers including 
both colour coherence and soft gluon interference effects. Hadronisation is simulated by the 
non-perturbative splitting of the gluons generated in the parton showers to create colour neutral 
clusters, followed by the fragmentation of those clusters. 



In addition to the above, DJANGO 6.0 [41| is used to study the effects of initial state pho- 



ton radiation. This Monte Carlo generator simulates the electroweak interaction using HERA- 



CLES [|2| and the hadronic final state using LEPTO. 



The DDIS Monte Carlo used is RAPGAP 2.1 |43]. In this generator a pomeron flux, de- 
pendent only on t and xjp, is associated with the proton. The pomeron is assumed to have 
partonic structure which evolves with according to the DGLAP prescription. The virtual 
boson interacts with a parton within the pomeron, the matrix element used being as for the ND- 
DIS Monte Carlo generators. Further parton emission is simulated using either parton showers 
(RAPGAP-PS) or the Colour Dipole Model (RAPGAP-CDM) and hadronisation is performed 
using JETSET. In conjunction with RAPGAP only the default fragmentation parameters are 
used. The pomeron parton distributions chosen give a good description of the latest HI mea- 



surements of F^^^^ 1 44 



Calculations of the expectations for the hadronic final state of instanton induced processes are 
made using QCDINS |^^, which is a version of HERWIG modified to generate QCD instanton 
induced events. The number of flavours generated in the instanton decay is determined by the 
available instanton sub-process centre-of-mass energy, y^. The decay is assumed to be isotropic 



in that frame. The minimum value of the so-called Holy Grail function [17|, which enters the 
expression for the cross-section for instanton induced events, is fixed to be the value it reaches 
at x' = 0.2, where x' = Q''^/{s' + Q'^) and Q'^ = —q'"^ is the negative invariant mass of the 
virtual quark in the instanton sub-process, see figure |l]. 
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Figure 1: Diagram illustrating the kinematics of instanton induced deep-inelastic positron-proton 
scattering. 

The number of Monte-Carlo events used in comparisons with the measurements presented 
in the following always exceeds the number of events in the data by a factor of three or more. 



4 Data Taking and Event Selection 

The data used here were collected in the 1994 running period in which 820 GeV protons and 
27.5 GeV positrons were brought into head-on collision in HERA. 

The trigger used in the collection of the data is the requirement that there be a localised 

energy deposit in the BEMC of at least 4 GcV with no associated veto signal from the ToF 
system. This has been shown to be more than 99% efficient for the collection of DIS events in 
which the scattered positron enters the BEMC and has energy greater than 10 GeV. 

The events are then required to contain a localised BEMC energy deposit of at least 12 GcV 
spatially associated with a signal in the BPC and with a shape compatible with being due to 
an electromagnetic shower. It is also demanded that at least one charged particle be detected 
in the CTD, allowing the determination of the z coordinate of the primary interaction vertex. 
This is required to be within 30 cm of the z coordinate of the expected vertex position. 

The polar angle of the positron candidate is required to lie in the range 156° < 9 < 173°. The 
range in the kinematic variables studied is restricted to W < < 70GeV^, 10~^ < a; < 10~^ 
and 0.05 < y < 0.6, all variables being determined using the measured polar angle and energy of 
the scattered positron. The lower restriction in the y range ensures that these determinations are 
of sufficient accuracy. The above criteria are such that the boundaries of the kinematic region 
studied are given by the restrictions imposed on the Lorentz invariant variables Q^, x and y. 
The number of events which satisfy the above is 53360. Less than 1% of these are background 
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events. The largest source of background is photoproduction interactions in which one of the 
final state hadrons fakes the required BEMC electron signal. 







The luminosity to which the data correspond is 1.32 pb , determined to a precision of 1.5% 
by studying the reaction ep — > ep'j using the luminosity measurement system. 

Similar criteria to those used to isolate DDIS events in the study of the structure of deep- 



inelastic diffraction |24] are used here to identify diffractive events. These are: the energy in 
the Plug must be less than 3GeV; there must be fewer than 3 hit pairs in the forward muon 
drift chambers and the most forward significant energy deposit detected in the LAr calorimeter 
must be at a pseudo-rapidity of 3.0 or less. Events in which the forward-going proton remnant 
is separated from the centrally produced hadronic system by a rapidity gap between rj ~ 6.6 
and T] = 3.2 and has a mass, uir, less than 1.6GeV/c^ satisfy these criteria. The DDIS event 
sample is restricted to the kinematic region defined by xjp < 0.05 and contains 3 491 events. 
The variable xp is calculated from the values of and determined using the parameters 
of the scattered electron and the value of Mx obtained from the measurement of the hadronic 
final state. 

As comparisons are made with Monte Carlo simulations which do not explicitly model the 
diffractive component of DIS, additional criteria are applied to select a sample of non-diffractive 
DIS events. These are that the events do not satisfy the diffractive selection above and in addition 
that the energy in the LAr calorimeter in the angular region 4.4° < 9 < 15° (3.26 > ij > 2.03) 
exceeds 0.5 GeV. The number of events in the NDDIS sample is 46 684. The non-diffractive 
Monte Carlo simulations are expected to model the hadronic final state of these events. To 
ensure that any events produced by these simulations which fall into the DDIS category are 
removed, before comparisons are made with NDDIS data the Monte Carlo events are required 
to have a hadronic energy of at least 0.5 GeV in the above angular region. 



5 Kg and A Identification 

Neutral kaons are identified through the decay 

i^^^vr+vr-, 

and A baryons through the decay 

A — > pTr~ 

and the charge conjugate reaction. The identification utilises the separation of the decay and 
primary vertices arising from the relatively long Kg and A lifetimes. 

The first stage in the identification of candidates for the above decays is a search for pairs 
of oppositely charged tracks in the CTD which originate from a vertex radially separated from 
the primary event vertex by at least 2 cm. Only well-measured tracks with pT > 0.15 GeV/c are 
considered. An attempt is made to fit these track pairs using the hypothesis that they originate 
from a common secondary vertex resulting from the decay of a neutral particle, generically 
termed a V^, which in turn originates from the primary event vertex. The candidates which 
satisfy this hypothesis are then required to have transverse momenta in the range 0.25 < < 
4.5(GeV/c)^ and pseudo-rapidities in the range —1.3 < t] < 1.3 to ensure that they are well 
within the acceptance of the CTD and that the mass resolution is good. 

In order to define a set of Kg candidates further requirements are made: 
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• Track pairs are required not to be A decay candidates. That is, it is required that the 
invariant mass of the track pair, mp-n^, satisfy m^^ > 1.125 GeV/c^, where the higher 
momentum track is assumed to be that of a proton and the lower that of a pion. 

• The difference of the distance of closest approach to the primary event vertex in the 
transverse plane. Ad, of the tracks making up the pair is required to satisfy 

1 + 10exp[-pT] 
|Ac!| > cm, 



where px (in GeV/c) is the transverse momentum of the candidate. The distance of 
closest approach of a track, d, is given a sign that depends on the track's charge and 
hence large values of |A(i| result for track pairs due to the decay of a into oppositely 
charged particles at a vertex removed from the primary vertex. The applied restriction is 
a function of the px to counteract the tendency of large px to cause small values of \d\. 

• It is required that | cos0*| < 0.95, where 6* is the angle in the transverse plane between 
the vector linking the primary and secondary vertices and the charged track direction in 
the rest frame of the track pair. 

A set of potential A decays is defined from the sample of candidates by demanding that: 

• Track pairs are not Kg decay candidates. That is, it is required that the invariant mass 
of the track pair, m^jr, satisfy m,r7r < 0.45 GeV/c^, assuming both tracks are due to pions. 

• Track pairs do not result from photon conversions. That is, the mass of the track pair 
rriee, assuming these to be electron and positron, is such that ruee > 0.05GeV/c^. 

• The distance of closest approach of the tracks making up the pair is required to satisfy 
|Ad| > 0.5 cm. 

The invariant mass spectra of the track pairs for the above samples are shown in figure 
These are calculated assuming that both tracks are due to pions for the K'^ candidates and for 
the A candidates that the higher momentum track is due to a proton, that with lower momentum 
to a pion. The final step in the selection of K'^ candidates is the requirement that the invariant 
mass of the track pair, niT^T^, satisfy 0.45 < m,r7r < 0.55GeV/c^. The A candidates are selected 
by requiring that mp-n^ satisfy 1.1066 < mpj^ < 1.1246 GeV/c^ . The numbers of Kg and A 
candidates are 1 813 and 395 respectively. 

The masses of the K^ and A are extracted from the spectra shown in figure || by fitting 
a Gaussian to the obvious peaks, using various functional forms to describe the remaining 
background. The backgrounds to the Kg and A signals are 8 ± 3% and 29 ± 7% respectively. 
The resulting masses are rrij^o = 497.9 ± 0.2 ± 0.2 MeV/c^ and niA = 1115.6 ± 0.4 ± 0.2 MeV/c^, 
where the first error is statistical and the second systematic. The systematic error reflects the 
effects of varying the assumptions on the functional form of the background remaining under the 
signal and of residual uncertainties in the calibration of the tracking system. These values are 
in agreement with the current world averages quoted by the Particle Data Group (PDG) p^ ], 
namely m^?^ = 497.672 ± 0.031 MeV/c^ and m^^^ = 1115.683 ± 0.006 MeV/c^. 

The numbers of Kg mesons and A baryons after background subtraction are 1 655 it 41 and 
237 lb 30 respectively. There are 116 it 21 A baryons and 126 it 20 A anti-baryons. These latter 
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Figure 2: Pion-pion (upper graph) and proton-pion (lower graph) mass spectra for candidates 
in the kinematic region 10 < < 70GeV^, 10""^ < x < 10~^ and 0.05 < y < 0.6 following the 
selection procedure described in the text. 
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numbers do not sum to give the total A plus A yield above as they are obtained from independent 
fits with independent background subtractions. 

The detection efficiencies for Kg mesons and A baryons in the given and r] ranges are 
determined using Monte Carlo techniques. The efficiency for detection is 25% at a pj^ of 
0.25 (GeV/c)^ and rises to 40% for a of 1 (GeV/c)^ and above. The detection efficiency is 
independent of pseudo-rapidity over the measured range. The A detection efficiency is 12% 
for all T] and py. There is no significant difference between the A and A detection efficiencies. 
Neither the Kg nor the A detection efficiencies depend significantly on x ov Q^. 

The Kg and A lifetimes obtained using the above event samples are r^o = 87.7 it 7.0 it 4.5 ps 
and TA = 330 ± 130 ± 30 ps respectively. The lifetimes are determined from the proper decay 
length distributions after subtracting the backgrounds extracted using the mass spectra. These 
are largest at small decay lengths. A correction is also applied for the decrease in K'^ and A 
detection efficiency which occurs at small decay lengths. This is determined using the LEPTO 
Monte Carlo with default JETSET parameters. The systematic error reflects the uncertainties 
in the background subtraction mentioned above and in the correction applied. The latter are 
estimated using half the spread of results obtained when the correction is derived using the 
LEPTO and CDM Monte Carlo programs with default JETSET parameters, with an additional 
contribution to cover the effects of deficiencies in the simulation of the HI apparatus. The 
good agreement of the above results with the PDG world averages, t^^"^ = 89.26 ± 0.12 ps and 

T^^^ = 263.2 lb 2.0 ps, confirms that the background subtraction and correction procedures are 
well understood. 



6 Results 

The observed numbers of Kg mesons and A baryons, after correction for event selection efficiency, 
particle detection efficiencies and branching ratios into the observed decay channels, correspond 
to the production of n^o = 16414 ± 403 ± 709 K^ mesons and ua = 2 879 ± 373 ± 341 A 
baryons in the pseudo-rapidity and transverse momentum range defined by —1.3 < ry < 1.3 
and 0.25 < < 4.5 (GeV/c)^ and in the kinematic region 10 < < JOGeV^, 10"^ < 
X < 10~^ and 0.05 < y < 0.6. The corrections are calculated by comparing the numbers of 
K^ mesons and A baryons in events generated using the LEPTO Monte Carlo with default 
JETSET parameters with the numbers found after a Monte Carlo simulation of the effects of 
the HI detector, reconstruction of the Monte Carlo data using the same programs as used for 
the HI data, and application of the selection criteria described above. The systematic errors 
represent half the spread in the determination of the correction factors arising from the use 
of the LEPTO and CDM generators with default JETSET parameters, the aforementioned 
uncertainties in the parameterisation of the background mass spectra in the Kg and A mass 
regions, residual uncertainties in the calibration of the tracking system and deficiencies in the 
simulation of the HI apparatus in the Monte Carlo program used. The effect of hard initial 
state photon radiation on these numbers, studied using DJANGO, is less than 2%. Hence no 
correction is applied, rather a contribution is added to the systematic error to cover the slight 
uncertainties introduced by neglecting this effect. 

Monte Carlo simulations indicate that the above numbers correspond to about 20% of the 

total number of K^ mesons and the same proportion of all A baryons produced in the x and 
region studied. The systematic uncertainties associated with the extrapolation from the 
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measured rj and pj^ region to all t] and pj, are large. Hence results quoted in this section are for 
the restricted region in which accurate measurement is possible. 

For the reasons discussed in the introduction, it is of interest to study the and A spectra 
in rj and pj^. The BFKL mechanism may cause excess strangeness production at large values 
of ry and p'^, whereas instanton induced events contain an excess of strange particles in the 
range ^ r] <, 3. The relevant measured spectra for production in NDDIS, in the kinematic 
region defined above, are shown in figure ^. This figure has been corrected for acceptance 
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Figure 3: Measured rj spectrum (upper graph) for production in the range 0.25 < p"^ < 
4.5 (GeV/c)^ and pj^ spectrum (lower graph) in the range —1.3 < r] < 1.3 in non-diffractive 
DIS for the kinematic region 10 < Q"^ < TOGeV^, 10""^ < x < 10"^ and 0.05 < y < 0.6 
compared with predictions from LEPTO with default JETSET parameters (LEPTO (Default)) 
and LEPTO with DELPHI JETSET parameters (LEPTO-D), see table 1. 
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and detection efficiency effects using Monte Carlo simulations, as have all subsequent figures. 
Here and in all following figures in which two sets of error bars are displayed the total errors 
are the result of adding in quadrature the statistical errors, shown using the inner error bars, 
and the systematic errors. The latter are estimated as previously described. Comparisons are 
made with the predictions of LEPTO using the default fragmentation parameters (A^ = 0.3), and 
using the set of parameters extracted by the DELPHI collaboration (As = 0.23, see table |l]). The 
normalisation of the Monte Carlo predictions is to the corrected number of NDDIS events in the 
specified x and range. In this and in other distributions it is observed that the description of 
the data provided using the DELPHI JETSET parameters in the LEPTO and CDM Monte Carlo 
models is significantly better than that obtained with the default settings. The parameter change 
that has the largest effect on the Monte Carlo predictions discussed here is the modification to 
As. A similar level of agreement is obtained if the strangeness suppression factor favoured by 
E665 (As = 0.2) is used. Once the modification to the fragmentation parameters has been 
made the description of the data is reasonable. There is no evidence for anomalous strangeness 
production in either the absolute normalisation or the shape of the measured distributions. 

The above agreement of the Monte Carlo predictions for the average number of mesons 
per event with the measurement does not prove that the distribution of the number of 
mesons per event is correctly described. The models might produce a larger number of strange 
particles in a smaller proportion of the events than observed in the data, or vice- versa. However, 
this is not the case, as may be seen from the following observation. The corrected proportion 
of NDDIS events containing two kaons in the kinematic region studied is 0.044 it 0.009 it 0.012. 
This number is to be compared with the LEPTO Monte Carlo predictions of 0.065 and 0.056 
with default and DELPHI JETSET parameters respectively, and with the corresponding CDM 
predictions of 0.060 and 0.044. HERWIG predicts that the proportion of di-kaon events is 0.059. 
Agreement with the Monte Carlo models which use JETSET is again observed to be better if 
the lower As value is used. There is no evidence for anomalous production of events containing 
many kaons. 

The A spectra in r] and are shown in figure ^ The statistical precision of these data does 
not allow discrimination between the default and DELPHI sets of JETSET parameters. The 
broad agreement with the Monte Carlo predictions indicates that the mechanisms responsible 
for di-quark production in the fragmentation in DIS at low x are similar to those at work in 
other processes and at higher values of x. 

In order to allow some separation between strangeness produced in the hard scattering and in 
the softer fragmentation process, production is also studied as a function of the fragmentation 
variable z, defined as 

_ + Vl,k^ 
IE, ' 

where i?^o and ^'^^ the energy of the and the component of its momentum along the 

struck quark direction, given by the four-momentum q + xp. This may be calculated using the 
measured parameters of the scattered positron and is the direction expected for the scattered 
quark in the naive Quark Parton Model. The quantity Eq is the energy expected for the struck 
quark, assumed massless, using the same approximation. Strange quarks may be produced in 
the hard scattering either through the interaction of the virtual boson with a strange sea quark, 
or directly via the BGF process, or through the decay of c quarks produced via BGF. The 
mesons formed from strange quarks produced by these means, referred to as "hard K^s^\ 
are likely to be at larger z than those produced in the fragmentation chain following, say, the 
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Figure 4: Measured rj spectrum (upper graph) for A production in the range 0.25 < Pt < 
4.5(GeV/c)^ and pj^ spectrum (lower graph) in the range —1.3 < rj < 1.3 in non-diffractive 
DIS for the kinematic region 10 < < TOGeV^, 10""^ < x < 10"^ and 0.05 < y < 0.6 
compared with predictions from LEPTO with default JETSET parameters (LEPTO (Default)) 
and LEPTO with DELPHI JETSET parameters (LEPTO-D), see table 1. 
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Figure 5: Measured z spectrum for production in the range 0.25 < v\ < 4.5(GeV/c)2 
and —1.3 < rj < 1.3 in non-diffractive DIS for the kinematic region 10 < < 70GeV^, 
10~^ < a; < 10~^ and 0.05 < y < 0.6 compared with predictions from LEPTO with DELPHI 
JETSET parameters (LEPTO-D) and the CDM with DELPHI JETSET parameters (CDM-D), 
see table 1; also shown is the spectrum obtained from LEPTO if only mesons which contain 
a strange quark produced in the hard interaction are considered (LEPTO (Hard K^)). 



interaction of the virtual photon with a u quark. Hence anomalous production at large z 
is likely to be associated with the hard sub-process, that at smaller z with the hadronisation. 
The K'^ spectrum in z is shown in figure |5|. The Monte Carlo curves displayed in this figure 
are obtained using LEPTO and the CDM. Of interest is the curve obtained when only hard 
mesons are considered. These are found to populate the high z region of the distribution, 
as expected, and form about 20% of the total number of mesons produced by the LEPTO 
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Figure 6: Mean numbers of K mesons per NDDIS event produced in the range 0.25 <pI < 
4.5(GeV/c)^ and —1.3 < rj < 1.3 as a function of x, with the restriction 0.05 < y < 0.6, in 
non-diffractive DIS in the regions indicated compared with predictions from LEPTO with 
DELPHI JETSET parameters (LEPTO-D) and the CDM with DELPHI JETSET parameters 
(CDM-D), see table 1. 

generator in the kinematic region studied when using the DELPHI JETSET parameters. About 
40% of these result from c quark decays. There is no evidence for anomalous production of 
strangeness associated with either the high or low z regions of the distribution. 

In figure ^ are shown the mean numbers of mesons observed per NDDIS event as a 
function of x in the regions 10 < < 15 GeV^, 15 < < 35 QeV^ and 35 < < 70 GeV^. 
There is no evidence for any systematic deviation as a function of either or x between 
the Monte Carlo predictions shown here and the data. As stated in the introduction, BFKL 
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dynamics may cause additional production at large values of rj and at small x. It has been 
argued that the CDM gluon radiation prescription models to some extent the effects of BFKL 



evolution on the hadronic state [47|. The similarity of the LEPTO and CDM predictions in 
figure ^ would then indicate that, in the rj region studied, there is no sensitivity to any putative 
additional strangeness production due to BFKL dynamics. A definitive statement must await a 
Monte Carlo implementation of the BFKL evolution prescription. 

An alternative way of examining the production data is to consider the ratio of the 
number of mesons produced to the number of primary charged tracks. This latter is defined 
to be the number of primary charged tracks in the transverse momentum and pseudo-rapidity 
region given by pT > 0.15 GeV/c, —1.3 < r/ < 1.3 and is corrected for detection efficiencies using 
Monte Carlo simulations. This means of presenting the data maintains sensitivity to effects 
that change the ratio of strange to non-strange quark production while reducing the effect of 
global multiplicity changes. Hence, for example, shifts in the strangeness suppression factor 
may be observed with reduced sensitivity to hadronisation variables which tend to change the 
overall track multiplicity. The production ratios, in the regions defined above, are shown 
in figure |^. Again, better agreement with Monte Carlo simulations is observed when using the 
DELPHI than when using the default JETSET parameters, providing further evidence that it 
is indeed the suppression of strangeness production in the hadronisation that is responsible for 
the differences, rather than a global discrepancy in the modelling of the multiplicity. 

Comparisons of the spectra in rj and pj^ for DDIS events with calculations using RAPGAP 
are shown in figure |8[ These data are corrected to the range xjp < 0.05 and to the x and 
region defined above using the results of the RAPGAP-PS Monte Carlo and a full simulation 
of the effects of HI detector. The RAPGAP-CDM predictions, made with the default JETSET 
parameters, are seen to lie a little above the measurements, as is the case for NDDIS. The 
agreement with RAPGAP-PS is somewhat better. 

The total rates of strange particle production in both DDIS and NDDIS are listed in table ^. 
They reveal that the level of strangeness production in DDIS is somewhat lower than that in 
NDDIS. This effect is probably due to the decreased overall level of particle production in DDIS 
in the r] region studied, a consequence of the lower hadronic mass of the fragmenting system 
{Mx < W), rather than a difference in the strangeness suppression factor. This argument is 
supported by the observation that the ratio of mesons to all charged particles in DDIS is 
consistent with that in NDDIS. Given that the overwhelming contribution to proton structure 
in the x range studied here is made by gluons, this observation is consistent with the electron- 
gluon scattering picture of DDIS. The consistency with the picture in which the scattering occurs 
between the positron and a flavour singlet pomeron is demonstrated by the agreement observed 
between the RAPGAP predictions and the DDIS measurements of the to charged track 
ratio. 

It is also apparent from table || that the difference in the two RAPGAP DDIS calculations 
of the multiplicity largely disappears when the ratio of that multiplicity to the primary 
charged particle multiplicity is considered. This occurs as the RAPGAP-CDM generator tends 
to produce a higher charged particle multiplicity than that seen with RAPGAP-PS. 



7 Comparison with Previous Measurements 



Figure |9| shows the rj spectra for and A production in DIS (NDDIS and DDIS) after correction 
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Figure 8: Measured rj spectrum (upper graph) for production in the range 0.25 < < 
4.5 (GcV/c)^ and pj^ spectrum (lower graph) in the range —1.3 < r) < 1.3 in diffractive DIS for 
the kinematic region 10 < < TOGeV^, 10"^ < x < 10"^, 0.05 < y < 0.6 and xp < 0.05 
compared with predictions from RAPGAP-PS and RAPGAP-CDM. 
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graph) in the kinematic region 10 < < 640 GcV^, 0.0003 < a; < 0.01 and y > 0.04; the 
results are for the range 0.5 < pt < 4.0GeV/c and the A results for 0.5 < pT < 3.5GeV/c, 
compared with the measurements of the ZEUS collaboration made in deep-inelastic electron- 
proton scattering at HERA. 
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Data-Set 


No. i^'^/event 


No. A/event 


No. K^/No. tracks 


HI NDDIS data 


0.278±0.008±0.012 


0.053±0.007±0.007 


0.058±0.002±0.003 


CDM (default) 


0.316 


0.052 


0.066 


CDM (DELPHI) 


0.284 


0.047 


0.062 


CDM (E665) 


0.260 


0.042 


0.054 


LEPTO (default) 


0.320 


0.055 


0.068 


LEPTO (DELPHI) 


0.282 


0.047 


0.064 


LEPTO (E665) 


0.267 


0.040 


0.056 


HERWIG 


0.316 


0.047 


0.068 


HI DDIS data 


0.238±0.026±0.016 


0.026±0.010±0.025 


0.063±0.007±0.004 


RAPGAP-PS 


0.255 


0.032 


0.072 


RAPGAP-CDM 


0.316 


0.040 


0.074 



Table 2: Measured mean and A multiplicities with 0.25 < < 4.5 (GeV/c)^ and -1.3 < 
rj < 1.3 for the kinematic region 

10 < Q2 < 70GeV^ 10"^ <x < 10"^ and 0.05 < y < 0.6 and 
the ratio of the above multiplicities to the mean multiplicity of primary charged particles with 
Pt > 0.15GeV/c in the same kinematic region in diffractive and non-diffractive DIS together 
with the predictions of various Monte Carlo models. 



to the transverse momentum range 0.5 < pT < 4.0GeV/c in the case of the mesons and 
0.5 < Pt < 3.5GeV/c for the A baryons and the kinematic region 10 < < 640 GeV'^, 
0.0003 < X < 0.04 and y > 0.04. The correction is performed using the LEPTO Monte Carlo 
with DELPHI JETSET parameters. The systematic errors shown include contributions due 
to the uncertainties in the correction procedure estimated using the spread of results obtained 
when the default JETSET parameters are used with LEPTO and both the default and DELPHI 
JETSET parameters with the CDM. Although the central values obtained are very similar to 
those shown in figures ^ and the correction is to a kinematic region which includes a significant 
range outside the reach of the current measurement and hence the associated systematic errors 
are larger than those shown in figures ^ and ^. These results are directly comparable with the 
measurements of the ZEUS collaboration also shown in figure ^, obtained in deep-inelastic 
electron-proton scattering at HERA. The two sets of data are seen to be in reasonable agreement. 

In order to facilitate comparison with other measurements, in the remainder of this section 
the and A NDDIS production results are presented as functions of various quantities in the 
hadronic centre-of-mass frame with the z direction chosen to be that of the exchanged virtual 
boson. The transformation to this frame is made using the measured energy and angle of the 
scattered electron. 



Figure IC shows the production rate as a function of Feynman-x, the fractional longitu- 
dinal centre-of-mass momentum xp = 2p'^/W. The data are also given in table ^. The xp range 
shown is restricted to the region in which, for the laboratory pseudo-rapidity and transverse 
momentum region over which particle identification is possible, the acceptance is at least 20%, 
namely 0.05 < xp < 0.40. The average acceptance over this range is 35%. The resulting mean 
centre-of-mass energy is (W) = 138 GeV. The measurements are corrected for the detector ac- 
ceptance, the region of laboratory pseudo-rapidity and transverse momentum in which particle 
identification is performed, and the kinematic region studied. The correction is performed using 
the LEPTO Monte Carlo with DELPHI JETSET parameters. The systematic errors represent 
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Xp 


1 dn 
N dxp 


0.05 - 0.10 


3.06 ± 0.23 ± 0.83 


0.10 - 0.15 


1.60 ±0.14 ±0.19 


0.15 - 0.20 


1.07 ±0.11 ±0.15 


0.20 - 0.25 


0.68 ±0.09 ±0.12 


0.25 - 0.30 


0.48 ±0.08 ±0.16 


0.30 - 0.40 


0.28 ± 0.05 ± 0.07 



Table 3: The spectrum as a function of Feynman-x (xp) at (W) = 138 GeV for the range 
10 < Q2 < TOGeV^ (see figure 10). 



the dispersion arising from using the default JETSET parameters with LEPTO and both the 
default and DELPHI JETSET parameters with the CDM in the correction procedure and the 
effects of the accuracy with which the transformation to the centre-of-mass frame can be per- 
formed. The shifts resulting from initial state photon radiation, studied using DJANGO, are 
found to be small and hence are not corrected for. The systematic error contains a contribu- 
tion covering their effects. Both the CDM and LEPTO Monte Carlo models describe the xp 
spectrum well, provided the DELPHI JETSET parameters are used. The agreement with the 
HERWIG Monte Carlo is also good. 



Also shown in figure M are the data of the EMC collaboration |48], obtained in muon-proton 



scattering at \/(W^ = 11.4 GeV, and of the E665 collaboration [^, obtained in muon-nucleon 
scattering at (W) = 17.1 GeV. These lie below the HI data and fall less steeply with xp, the 
differences being more pronounced for the lower energy EMC dataQ. This observation may be 
largely explained by the extra phase space available for QCD radiation at the higher centre-of- 
mass energy available at HERA. 

Comparisons of the xp spectra with electron-positron data are complicated by the direct 
production of heavy quarks in electron-positron annihilation, with subsequent decay to strange 
quarks and hence enhanced strange particle production. This problem is partially solved here 
by making the comparison with Monte Carlo predictions for electron-positron annihilation to u, 
d and s quarks only, at the HI centre-of-mass energy. Residual differences are to be expected as 
the proportion of primary s quarks produced in the simulation is still not the same as that in 
DIS at this energy. The Monte Carlo used (JETSET with DELPHI tuning) accurately describes 



electron-positron annihilation data when the production of all possible flavours is allowed |11]. 
A further complication is that the electron-positron annihilation data are available in terms of 
the fractional centre-of-mass momentum, Xp = 2p/W. However, the transverse momenta in the 
hadronic centre-of-mass frame, p^, of the mesons studied here are generally much smaller 
than their longitudinal momenta and so the momentum p ^ pp. Hence xp = Xp io sea accuracy 
better than that represented by the symbol size in the figures. The spectrum in Xp from 
the JETSET simulation, divided by two to compensate for the production of a quark anti-quark 
pair in electron-positron annihilation, may thus be expected to be similar to the xp spectrum 
resulting from the struck quark in NDDIS. That this is the case is seen in figure 

Despite their lower precision, the A measurements, shown in table |^ and figure 0, reveal a 



^Note that the EMC data here are at a lower centre-of-mass energy than those compared with HI results 
in il. 
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Figure 10: Production of mesons as a function of Feynman-x {xp) in the range 10 < < 

70 GeV'^ compared with the muon-proton and muon-nucleon scattering results of the EMC and 
E665 collaborations, the results of a Monte-Carlo simulation of electron-positron annihilation to 
u, d and s quarks at the same centre-of-mass energy as the HI data (JETSET PS) and various 
other Monte Carlo simulations of DIS. 
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Figure 11: Production of A baryons as a function of Feynman-a; (xp) in the range 10 < < 

70 GcV^ compared with the muon-proton and muon-nucleon scattering results of the EMC and 
E665 cohaborations, the resuhs of a Monte-Carlo simulation of electron-positron annihilation to 
u, d and s quarks at the same centre-of-mass energy as the HI data (JETSET PS) and various 
other Monte Carlo simulations of DIS. 
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Xp 


1 dn 
N dxp 


0.05 - 0.15 


0.46 ±0.10 ±0.10 


0.15 - 0.25 


0.13 ±0.05 ±0.03 


0.25 - 0.45 


0.12 ±0.03 ±0.04 



Table 4: The A spectrum as a function of Feynman-x at {W) = 138 GeV for the range 10 < 
Q2 < 70GeV2(see figure 11). 



XF 




0.05 - 0.10 


0.47 ±0.06 ±0.09 


0.10 - 0.15 


0.84 ±0.06 ±0.38 


0.15 - 0.20 


0.85 ±0.09 ±0.27 


0.20 - 0.40 


1.70 ±0.08 ±0.42 



Table 5: Corrected values of {pr) ^ ^ function of Feynman-x at (VF) = 117 GeV for the 
range 10 < < roGeV^ (see figure 12). 



similar trend to the K results discussed above when compared with the lower energy EMC 1 45 ] 
and E665 measurements. The agreement with the results of a simulation of electron-positron 
annihilation as described above is good. The HERWIG expectations differ from those of LEPTO 
and the CDM somewhat, but the data are not precise enough to discriminate between the models. 

The relationship between the fractional longitudinal momentum and the mean squared trans- 
verse momentum, which is sensitive to QCD radiation effects, is often illustrated in the "seagull" 
plot. Such a plot is shown for production in figure |l^, the data also being given in table ^. 
Here, the W range of the data is restricted in order to reduce the mean centre-of-mass energy 
to (yV) = 117 GeV, comparable with previously published HI data on charged particle produc- 
tion [^9|, also shown in the figure. The and charged particle measurements are seen to be 
in good agreement and both lie at higher (p^?) for a given xp than the lower centre-of-mass 
energy measurements of the E665 and EMC collaborations. This difference is expected within 
the framework of QCD due to the effects of increased radiation at the higher energy. The CDM 
and LEPTO models, with DELPHI JETSET parameters, are able to describe the data well. 
HERWIG (not shown) describes the data shown here with a similar level of precision. 

In order to study more directly the behaviour of and A production with the hadronic 
centre-of-mass energy, the mean and A multiplicities are shown as a function of W"^ in 
figure ^ and are also given in tables |^ and ^. These multiplicities are corrected so that 
they represent the numbers of particles produced with Xi? > 0. A procedure similar to that 
described above is used to obtain the correction and the associated systematic error. Both the 

and A measurements show that the logarithmic increase in mean multiplicity observed at 
lower energies (here data are shown from the E665 collaboration and from the anti-neutrino 
and neutrino-proton scattering experiments of the WA21 collaboration |8|) persists to the much 
higher energies available at HERA. 



27 



1 

0.9 
0.8 
0.7 
0.6 
0.5 

0.4 
0.3 

0.2 



0.1 



• H1 (K°), <W>=117GeV 

* H1 (charged particles), <W>=117 GeV 

❖ E665 (K°), <W>=17.1 GeV 
□ EMC(K*"), <WV^^=11.4 

CDM-D 
LEPTO-D 




0.05 



0.1 



0.15 



0.2 



0.25 



0.3 



0.35 



0.4 



0.45 



Figure 12: The mean squared transverse hadronic centre-of-mass momentum {{px)) of K 
mesons and of all charged particles as a function of Feynman-x {xp) in the range 10 < < 
70 GcV^ compared with measurements of and production in muon-nucleon and muon- 
proton scattering made by the E665 and EMC collaborations; also shown are various Monte 
Carlo predictions. 
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Figure 13: The mean (upper graph) and A (lower graph) multipUcities for Xi? > as a 
function of in the range 10 < < 70 GeV^ compared with the muon-nucleon scattering 
results of the E665 collaboration and WA21 measurements from anti-neutrino and neutrino- 
proton scattering. 
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W (GeV^) 




6400 - 19300 


0.74 ±0.03 ±0.09 


19300 - 32200 


0.93 ±0.05 ±0.11 


32200 - 45100 


0.86 ±0.06 ±0.08 


45100 - 58000 


0.82 ±0.09 ±0.16 



Table 6: Mean multiplicity for rcj? > in the hadronic centre-of-mass frame as a function of 
W"^ (see figure 13). 



W (GeV^) 


(A) 


6400 - 15625 


0.13 ±0.03 ±0.03 


15625 - 58000 


0.17 ±0.03 ±0.04 



Table 7: Mean A multiplicity for xp > in the hadronic centre-of-mass system as function of 
W"^ for the range 10 < < 70GeV2(see figure 13). 



8 QCD Instanton Production 

If a significant proportion of DIS events were induced by QCD instantons, large changes in 
the strangeness composition of the hadronic final state would be expected. This is illustrated in 



figure 14 where, to ensure reliable simulation using the instanton generator, the kinematic region 
studied is further restricted to 10 < < roQeV^, 10"^ < x < 10"^ and 0.1 < y < 0.6. The 
number of events in this region is 16 486 after the application of the NDDIS selection criteria. 
In order to determine an upper limit for the cross-section for instanton induced events it is then 
assumed that the measured rate of production, shown as a function of rj in figure 15, is 
the result of a proportion / of instanton induced events with the rest being due to the NDDIS 
models discussed above. A minimisation procedure is used to determine the value of /. 

When using LEPTO and the CDM to calculate the rate expected for NDDIS without an 
instanton contribution, the DELPHI value and other fragmentation parameters are used as 
these are determined in electron-positron annihilation, in which there is no contribution from 
instanton induced events. The results obtained are / = —0.044 ± 0.032 and / = —0.002 ± 0.031 
using QCDINS with HERWIG and the CDM respectively. The largest value of / is obtained 
using QCDINS with LEPTO and is / = 0.006 ± 0.030. The fit to the measured K° production 



rate resulting in this latter case is shown in figure 15. The / values are consistent with being 
zero and hence a 95% confidence level upper limit of 0.9 nb is placed on the cross-section for 
instanton induced events in the above kinematic region, where the largest / value is used in 
the determination of the limit. The limit is of the same order of magnitude as a recent, albeit 
uncertain, determination of the cross-section for QCD instanton induced events in DIS [50|. 



9 Conclusions 



Strangeness production in deep-inelastic positron-proton scattering (DIS) is studied using the HI 
apparatus at the HERA collider at DESY. The rates of meson and A baryon production are 
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Measured region 




Figure 14: Numbers of mesons per event with 0.25 < pj^ < 4.5 (GeV/c)^ as a function of rj 
for the kinematic region 

10 < Q2 < TOGeV^, 10-3 < x < 10"^ and 0.1 < y < 0.6 predicted by 
models of non-diffractive DIS (broken hues) and by QCDINS (continuous hne). 
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Figure 15: Measured numbers of mesons per event with 0.25 < pj^ < 4.5(GeV/c)^ for the 
kinematic region 10 < < 70GeV^ 10"^ <x < lO'^ and 0.1 < y < 0.6 (points) with the fit 
result (dotted line) and the fraction of instanton induced events (solid line), / = 0.006, obtained 
as described in the text; the errors shown are the statistical and systematic errors added in 
quadrature. 
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measured in the kinematic region 10 < < 70GeV^, 10~^ < x < 10~^ and 0.05 < y < 0.6. No 
evidence for anomalous sources of strangeness is observed. A reduction of the strangeness sup- 
pression factor, A^, from the current default value in the JETSET string hadronisation scheme, 
0.3, to a value of 0.23 or 0.2 significantly improves the description of the measurements when 
using both the LEPTO and CDM Monte Carlo models. This observation is in agreement with 
the results of previous measurements made in DIS |^ and electron-positron annihilation |11]. 



The agreement with the latter suggests that the colour forces involved in the hadronisation 
process are similar in electron-positron annihilation and DIS in the low x region studied here. 
Further evidence for this is the reasonable description of the measured A production rate given 
by Monte Carlo models using JETSET hadronisation. Such agreement has also been observed 



in electron-positron annihilation [51|, suggesting that the di-quark production mechanisms are 
similar in both processes. 

No evidence is found for differences in the relative rate of strangeness production in diffractive 
as opposed to non-diffractive DIS, an observation consistent with both the hypothesis that these 
events result from deep-inelastic positron-pomeron scattering and the hypothesis that they are 
the product of boson-gluon fusion reactions with subsequent modification of the quark and 
anti-quark colour charges such that these form a colour singlet. 

Direct comparisons of the fractional longitudinal hadronic centre-of-mass momentum [xp = 
2p*^/W) spectra with previous studies of and A production in DIS |48, ^ at lower energy 



reveal an increase in the number of particles, particularly at low function of energy. This 

is expected within the framework of QCD and may be explained as being due to the increased 
amount of radiation at the higher energy. The xp spectra are also compared with Xp = 2p/W 
spectra obtained from a simulation of electron-positron annihilation to u, d and s quarks at the 
same centre-of-mass energy. These are in good agreement, providing further evidence that the 
processes involved in the hadronisation in DIS at the low values of x accessible at HERA are 
the same as those in electron-positron annihilation. A further consequence of the increase in 
QCD radiation with energy is observed when the mean squared transverse momentum of the 
i^" mesons in the hadronic centre-of-mass frame {{pr)) is studied as a function oi xp. The 
high energy data are seen to be at higher {p^) for a given xp than the lower energy data |^, ^]. 
These observations are described well by models of DIS which incorporate QCD effects. 

As no evidence for anomalous levels of strangeness production is found, an upper limit of 
0.9 nb at the 95% confidence level is placed on the cross-section for QCD instanton induced 
events in DIS. 
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